Kinetic investigation on extrinsic spin Hall effect induced by skew scattering 



00 
O 

o 

(N 
C 



J. L. Cheng and M. W. Wtfl 

Hefei National Laboratory for Physical Sciences at Microscale and Department of Physics, 
University of Science and Technology of China, Hefei, Anhui, 230026, China 
(Dated: February 5, 2008) 

The kinetics of the extrinsic spin Hall conductivity induced by the skew scattering is performed 
from the fully microscopic kinetic spin Bloch equation approach in (001) GaAs symmetric quantum 
well. In the steady state, the extrinsic spin Hall current /conductivity vanishes for the linear-k 
dependent spin-orbit coupling and is very small for the cubic-k dependent spin-orbit coupling. The 
spin precession induced by the Dresselhaus/Rashba spin-orbit coupling plays a very important role 
in the vanishment of the extrinsic spin Hall conductivity in the steady state. An in-plane spin 
polarization is induced by the skew scattering, with the help of the spin-orbit coupling. This spin 
polarization is very different from the current-induced spin polarization. 

PACS numbers: 72.25.Pn, 72.25.Rb, 71.70.Ej, 71.10.-w 
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Generating and manipulating the spin polarization in 
semiconductors is one of the important prerequisites for 
the realization of the new spintronic device^ Spin Hall 
effect (SHE) is considered as the convenient method to 
generate spin polarization in additional to the traditional 
methods such as the external magnetic field, the cir- 
cular/linear polarized laser ^ the spin-galvanic effect^ 
and the spin injection from the ferromagnetic metal to 
semiconductor £ The SHE is induced by intrinsic or ex- 
trinsic spin-orbit coupling (SOC) 6 which gives rise to the 
spin current vertical to the charge current without ap- 
plying external magnetic field and/or spin accumulation 
at sample edgesi 7 -^ Experimentally, the spin Hall con- 
ductivity (SHC) is estimated indirectly by the spin accu- 
mulation at the sample edges ; 10 ' 11 : 12 ' 13 : 14 or the charge 
current with a transverse magnetic field applied in a 



gyrotropic systems 



15.16.17 



Recently, a direct electronic 



measurements of the SHE is given by Valenzuela and Tin- 
kham in the metallic conductor. 18 All these effects arc 
explained as the intrinsio^^^^i^^^^ 5 .^^ and/or 
the extrinsic SHE&L2£^&2L22£2£i&2S£L2& theoreti- 
cally by using the Kubo formul a 19 : 20 ' 30 : 31 or the Boltz- 
mann equatio n 6 : 7 : 23 : 34 : 36 with only the carrier-impurity 
scattering included. 

The intrinsic SHE is induced by the intrinsic SOC 
(i.e., the Dresselhaus^ 9 . and/or the Rashba 40 SOCs) with 
the applied external electric field and the resulting SHC 
was at first thought as dissipationless one in perfect 
crystals! 19 i 20 Later investigations proved that this kind 
of SHC disappears even for the infinitesimal impurity 
density with the vertex correctio n 22 : 23 : 41 : 42 : 43 : 44 in the 
Rashba model or the linear Dresselhaus model in quan- 
tum wellsj 4 ^ but remains a finite value for the cubic 
SOC. 46 However, the spin current is not an observable 
quantity. The Laughlin's gauge gedanken experiment in- 
dicates that the intrinsic SHE cannot lead to any spin ac- 
cumulation at sample edges ; 24 : 25 unless in a mesoscopic 
system. In additional to the Dresselhaus and the Rashba 
SOC, the mixing between the valance and the conduc- 
tion bands gives rise to two corrections: One is the 



additional spin-dependent electron-impurity or the 
electron-phono n 16 : 17 skew scattering. The extrinsic SHE 
induced by the skew scattering alone has been widely 
studied by using both the Kubo formula and the Boltz- 
mann equation method; 6 -^ : 28 : 29 : 30 : 34 : 35 ' 36 and the nonzero 
extrinsic SHC is obtained. Lately, Tse and Das Sarma 31 
have proved the vanishment of the extrinsic SHC by 
considering the vortex correction of the linear SOC in 
the Kubo formula. However, a fully microscopic cal- 
culation of the extrinsic SHC from the kinetic equation 
approach is still missing and the vanishment of the ex- 
trinsic SHC by the vortex correction from the Kubo ap- 
proach needs to be verified from the kinetic approach. 
The other is the additional spin-dependent position and 
velocity operators which bring the correction to the defi- 
nition of the spin current, and are referred to as the side- 
jump mechanisnifS. Here two corrections on the definition 
should be specified. The first comes from the electrical 
potential which gives an intrinsic-like contribution and 
again cannot contribute to the spin accumulation accord- 
ing to Refs. [2412511471 . The second comes from the spin 
dependent scattering in which high order correlations be- 
tween different wave vectors need to be considered It 
is hard to include the second correction in the Boltzmann 
equation approach, 6 though it also gives rise to the ex- 
trinsic SHE, together with the skew scattering. In the 
following, we concentrate on the extrinsic SHE induced 
by the skew scattering. 

The spin polarization is not only accumulated at the 
sample edges due to the extrinsic SHE, but also observed 
simultaneously inside the samples which is induced by 
the charge current j 10 ' 11 : 13 ' 49 : 50 ' 51 : 52 By comparing the ex- 
periments in Refs. [T^fl3ll5ll . it is easy to find that both 
spin polarizations are in the same order. In theory, En- 
gel et a/4 9 . and Trushin and Schliemann 5 ^ attributed the 
current-induced spin polarization (CISP) in the homoge- 
neous system as the results of the current induced effec- 
tive magnetic field (EMF) from the SOC— Tarasenko 
showed that the spin-flip phonon scattering in asym- 
metric quantum wells can also induce this kind of spin 
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polarization. 50 However, the spin-conserving skew scat- 
tering in two dimensional GaAs semiconductor can also 
induce spin polarization inside the sample due to the ex- 
trinsic SHE. This effect has not been studied in the liter- 
ature. Moreover, a fully microscopic kinetic investigation 
on the extrinsic SHE is also missing in the literature. 

In this paper, we focus on the kinetic process of the 
extrinsic SHE induced by the skew scattering in symmet- 
ric GaAs (001) quantum well from the kinetic spin Bloch 
equation (KSBE) approach) 53 ' 54 We demonstrate the im- 
portant role of the spin precession induced by the (intrin- 
sic) Dresselhaus/Rashba SOC to the SHE and show that 
it is inadequate to study the extrinsic SHE from the Kubo 
formalism without considering the Dresselhaus/Rashba 
SOC in the literature. We further show that the extrin- 
sic SHE can generate spin polarizations in homogeneous 
system. 

By using the non-equilibrium Green function method 
and the generalized Kadanoff-Baym Ansatz^ 5 - we con- 
struct the KSBE 5 ^ 4 - for electrons as follows 

dt 



dt dk x 
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Here, p k {t) 



is electron density ma- 
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trix with wave vector k at time t. The applied elec- 
tric field E is assumed along the cc-axis and the mag- 
netic field B is along the x-y (well) plane. The coher- 
ent term describes the spin precession along the mag- 



netic/effective magnetic field and is given by 
i[^(n D (k)+gp B B) -a, Pk (t)}, where n D (k)=^(k x (k 



coh 



(a ) 2 )>M(^) 2 - fe x)>°) represents the EMF from the 
Dresselhaus SOC22, with 7 standing for the material- 
determined SOC strengt h 56 i 57 and a being the well width. 
Infinitc-wcll-dcpth assumption is adopted here and only 
the lowest subband is taken into account due to small 
well width. The spin conserving scattering ^ 1 is 

scat 

given by the spin conserving electron-impurity scattering, 
the electron-phonon scattering and the electron-electron 
Coulomb scattering which are given in detail in Ref. l58l . 

The skew scattering is given by the third order expan- 
sion of the electron-impurity scattering and reads £ i 30 i 31 



dpi, 
dt 



kik 2 ;qi 92 

xV(k-k 1 ,q 1 )V{k 1 -k 2 ,q 2 )V{k 2 - 
x{(k - ki) x (k 2 - k) • <r,p k2 } . 



~ £k)<5(Ek! - £ k) 

k, -q 1 - q 2 ) 

(2) 



Here A7 = 9m l E (3-,,) is the strength of the skew scat- 
tering with 77 = A _rg £ Eg and A are the band gap and 
the spin-orbit splitting of the valance band separately. 

v(d, qz ) = ^fi^im) with i{x) 



_ e inx sin(7ra) 

7TX(1 — X 2 ) 



l 2 +q'j+K 2 

being the form factor and k denoting the Debye-Hiicke 
screening constant. This term produces an asymmetric 
spin-conserving scattering of electrons so that the spin- 
up/down electrons prefer to be scattered to the left /right 




FIG. 1: (Color online) The schematic for the skew scatter- 
ing and the intrinsic spin orbit coupling in the k-space. 
((g)) stands for the spin-up (-down) polarization. The blue ar- 
rows give the direction of the EMF from the Dresselhaus SOC 
which rotates the spin polarization to the direction indicated 
by the red dashed arrows. 



side. A schematic for the skew scattering is given in Fig. 
[T]for a left moving electron. 

We first analysis the KSBE for some simple cases. Af- 
ter carrying out the summation over k, one gets the 
equation of continuity for the spin density S = ^ k Sk = 
X. k Tr[/9 k cr] as 



8S(t) 
dt 



-^2fl D (k) x s k (t) -gp B B x S(t) = . (3) 



The summation of the spin-conserving scattering is natu- 
rally zero. The summation of the skew-scattering in Eq. 
([2]) for the two-dimension system is also zero because it 
can only skew the spin-up and -down electrons separately 
instead of flip them. Equation ([3]) is consistent with the 
result in Ref. |4J which is obtained by using the general 
operator commutation relations. When only the linear- 
k terms in f2 D (k) is retained, for the steady state one 
obtains 



J?, = 



eg (i B By 



y 7 m*((7r/a) 2 - <fc 2 » 



S z 



(4) 



Here J z = — e^ k ftk/m*(/kj — /ki) is the spin cur- 
rent which is the simplest but the most widely-used 
definition 6 without considering contribution from the 
side-jump effect. The SHC is given as er* = Jy/E. In the 
derivation, we also take the approximation ^ k k x kyS^ — 

-^{kl)J> with (kl /y ), 



'(k y )Jx an d J2kkykx s \<L 



the average value of ki/.. If the strain-induced SOC 

x/y 

H s = a(k y a x — k x tjy) (with the same expression of the 
Rashba SOC) is taken into account, Eq. (TJJ changes into 



J: 



l{{l?-{kl))B y -2aB x 



egp-B 



S z . (5) 



7 2 ((f) 2 -<fc 2 ))((f) 2 -(fc 2 ))-4a 2 m* 

Without the external magnetic field, Eqs. (H]) and |(SJ| give 
Jy = 0, which verifies the zero extrinsic SHC given by 
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Ref. [3l[ from the Kubo approach. However, the skew 
scattering does generate the spin currents when the elec- 
tric field is applied (in Fig. [I]), but it is just the dynamic 
one. To make clear how the spin current disappears, we 
rewrite the x-component of Eq. ([3]) as 



10 l 



dS x /dt = m* 7 [(^/a) 2 - (k 2 z )} J*/e 



(6) 



It is hence easy to find that the spin-Hall current is con- 
verted to the spin polarization along the x-axis and it 
tends to zero in the steady state. This can further be 
understood from Fig. [TJ after the spin current is excited 
by the skew scattering, the spin polarization distributes 
(0 and ®) anti-symmetrically at ±k y . However the y 
components of the EMF due to the SOC (the blue ar- 
rows) have the same symmetry and rotate the spin po- 
larizations to the ~x direction (the red dashed arrows). 
Therefore, the spin Hall current is converted to the in- 
planc spin polarization. 

Now we show the time evolution of the SHC afj and 
the spin polarization P x — S x /N e calculated by numeri- 
cally solving the KSBE with all the scattering explicitly 
included at T — 200 K in Fig. The parameters in the 
calculation are taken as following: the electron and im- 
purity densities N e = N% = 4 x 10 11 cm -2 ; a = 7.5 nm; 
E = 0.1 kV/cm; 7 = 11.4 eV-A 3 ; and y/X'j = 2.07 A. 
From the figure, the SHC increases with time first from 
zero value to a maximum one in nearly 1 ps, then de- 
creases slowly to a very small value (in stead of zero due 
to the inclusion of the Dresselhaus SOC with the cubic 
k terms) in a characteristic time scale about 50 ps. The 
spin polarization is along the — x-axis and increases from 
zero to its steady state P x — 1.2 x 10 -4 . 

The above evolution is easy to understand with the 
help of the schematic in Fig. [TJ When the positive elec- 
tric field is applied, the skew scattering scatters the spin- 
down electrons to the upper panel of the k-space and the 
spin-up ones to the lower panel. This leads to the spin 
currents flowing along the y-direction. The strength of 
the skew scattering is determined by the shift of the elec- 
tron distribution, so the SHC increases fast to its max- 
imum value o z y ~ 2.5 x 10 at the time scale of the 
charge current (see the dashed curve in the inset of Fig. 
[2]). Then the spin polarization precesses along the y- 
components of the EMF from the SOC to the —x direc- 
tion to generate the in-plane spin polarization. Due to 
the symmetry, the spin polarizations along the y and z 
axes are zero. 

It is interesting to see that the obtained SHC in the 
steady state is orders of magnitude smaller than the one 
obtained from the Kubo formula widely used in the litera- 
ture where only the lowest order diagram is considered. 30 
With only the impurity scattering included at zero tem- 
perature, the Kube formula gives the extrinsic SHC 3 ^ as 
2n m '\y 2 s F 0x _ 3 37x iQ-3 ax with thc charge conductiv- 
2 

ity cr x = This value is orders of magnitude larger 

than our result u z y ~ 4.1 x lO -8 ^, obtained from the fully 
microscopic KSBE. The difference is caused by the inho- 
mogeneous broadening 59 in spin precession due to the 
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FIG. 2: Time evolution of the SHC with (solid curve) /without 
(chain curve) coherent terms and the time evolution of the 
spin polarization (dashed curve) along the z-axis at T = 200 
K. TV, = N e = 4 x 10 11 cm -2 and the electric field E = 
0.1 kV/cm. The corresponding time evolutions of the SHC 
of the first 4 ps are shown in the inset, together with the 
charge conductivity o x . It is noted that the scales of the spin 
polarization/charge conductivity are on the right hand side 
of the frames. 



Dresselhaus/Rashba SOC. To show this effect, we drop 
the coherent term and plot the time evolution of thc SHC 
by chain curves in Fig. O We find that the steady SHC 
a* ~ 0.4 x lO -3 ^ is much closer to the one above given 
by the Kubo formula at K. Therefore, we conclude that 
it is not suitable to calculate the extrinsic SHC without 
considering the spin precession. Furthermore, the side 
jump mechanism gives an time-independent contribution 



as cjy = -2X^ 2 e 2 N e ~ -2 x 10 3 e 2 /h with the opposite 



sign of the one from the skew scattering. 

It is further noted that although the spin precession in- 
duced by the EMF leads to the vanishment of the SHC, it 
plays a very important role in generating the spin polar- 
ization. From the discussion above, it is easy to see that 
the skew scattering alone cannot generate any spin polar- 
ization due to the anti-symmetrical spin polarization at 
±k y . Differing from the CISP— where the SOC is used to 
provide a current induced EMF ) 10 i 58 here the SOC acts 
as an anti-symmetrical precession field, which is the same 
as the spin polarization induced by the spin-dependent 
phonon scattering. 3 We further note that the spin polar- 
ization induced by the skew scattering is very different 
from the CISP: (i) The spin polarization induced by cur- 
rent induced EMF prefers the small well width which 
gives a large EMFji 9 . whereas the one induced by the 
skew scattering and the spin precession prefers a large 
well width. Thc later can be seen as following: from 
Fig. O the evolution of the extrinsic spin Hall concurrent 
can be written as Jy(t) = JyQe~ t ^ TE with the relaxation 
time t s . J z q approximates to the maximum value which 
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is only determined by the skew scattering. Therefore 
S x = m*T s j{{n/a) 2 - (^))J* >0 /e. For the D'yakonov- 
Perel' mechanism^ r s oc [7((7r/a) 2 — (fc 2 ))]~ 2 - Hence 
S x cx 7((7r/a) 2 - (fc 2 ))" 1 . (ii) When the electric field 
is along the x-axis, for the Dresselhaus SOC, the CISP 
is along the x direction, while the polarization induced 
by the skew scattering is along the — x direction. How- 
ever, for the strain-induced or the Rashba SOC, the spin 
polarizations from both mechanisms are along the same 
direction, (hi) The CISP decreases with the impurity 
density^ because high impurity density reduces the EMF 
effectively However, the skew-scattering-induced spin 
polarization increases with impurity density as the skew 
scattering is proportional to the impurity density. 

In summary, we investigate the SHC and the spin po- 
larization induced by the k-asymmetric spin-conserving 
skew scattering in symmetrical (001) GaAs quantum 
well from the fully microscopic KSBE approach at high 



temperature with all the scattering explicitly included. 
We find the spin precession induced by the Dressel- 
haus/Rashba SOC has a very important effect on the 
extrinsic SHC and verify the vanishment of the SHC for 
linear k-dependent SOC. We also show that the SHC in- 
duced by the skew scattering calculated from the Kubo 
formula in the literature is inadequate without consider- 
ing the spin precession. Finally we show that with the 
joint effects from the skew scattering and the spin pre- 
cession, an in-plane spin polarization can be generated 
which can be further rotated to the z-direction by apply- 
ing an external in-plane magnetic field. 
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